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ELECTROSTATIC  ION-CYCLOTRON  WAVES  IN  MAGNETO  SPHERIC  PLASMAS: 

NON-LOCAL  ASPECTS 


I.  INTRODUCTION 

The  current  driven  ion-cyclotron  instability  (CDICI)  has  been  of 
considerable  interest  to  space  plasma  physicists.  Recent  observations  by 
Kintner  (1980)  and  Kintner  et  al.  (1978)  have  made  it  very  topical. 
Drummond  and  Rosenbluth  (1962)  were  first  to  examine  this  instability 
analytically  while  Kindel  and  Kennel  (1971)  studied  it  in  the  context  of 
space  plasmas  in  the  earth's  magnetosphere.  The  results  from  Kindel  and 
Kennel  have  frequently  been  applied  to  observations  involving  the 
electrostatic  CDICI  in  space.  It  should  be  noted,  however,  that  the 
analysis  of  both  Drummond  and  Rosenbluth  and  Kindel  and  Kennel  is  a  local 
analysis.  They  consider  a  uniform  zeroth  order  magnetic  field  (i.e.,  B_  = 
B02),  thereby  neglecting  the  self-consistent  magnetic  field  generated  by 
the  field  aligned  currents.  This  magnetic  field  (usually  small)  will  give 
rise  to  a  shear  in  the  zero  order  magnetic  field  (see  Figure  1)  and 
consequently  make  the  zero  order  field  space  dependent  (i.e.,  B(x)  = 

a  a 

B£x)y  +  B^z,  where  shear  is  in  the  x-direction) . 

In  general,  magnetic  shear  is  a  damping  agent  and  can  significantly 
alter  the  local  mode  structure  (e.g.,  it  was  shown  by  Ross  and  Mahajan, 
(1978)  that  an  infinitesimal  shear  can  completely  damp  the  universal  drift 
instability).  Recently  Waltz  and  Dominguez  (1981),  motivated  by  the  TFR 
(1978)  tokamak  experiment,  have  provided  numerical  results  for  the  behavior 
of  CDICI  in  a  sheared  magnetic  field  pertinent  to  the  TFR  (1978) 
parameters.  Ganguli  and  Bakshi  (1982)  have  given  a  detailed  analytical 
treatment  of  the  CDICI  in  a  sheared  magnetic  field  and  have  concluded  that 
even  a  small  shear  can  give  a  significant  reduction  in  the  growth  rates. 
It  also  greatly  reduces  the  region  of  unstable  perpendicular  wavelengths. 

The  above  mentioned  treatments  of  shear  implicitly  assume  that  a 
uniform  shear  prevails  over  a  large  enough  region  of  space.  A  finite 
current  channel  size,  of  width  Lc,  may  introduce  important  modifications 
when  Lc  becomes  comparable  to  or  smaller  than  the  shear  length  Lg. 
Recently  Bakshi,  Ganguli  and  Palmadesso  (1983)  have  developed  the  theory 
for  the  effects  of  shear  on  CDICI,  taking  into  consideration  the  role 
played  by  a  finite  width  current  channel.  The  main  result  of  that 
Manuscript  approved  August  9,  1983. 


investigation  shows  that  the  growth  rate  is  given  by  the  pure  shear  theory 
when  Lc  >  Lg  and  reverts  to  the  local  growth  rate  if  Lc  <<  Lg,  with  a 

smooth  transition  as  a  function  of  (Lc/Lg).  There  is  another  non-local 
effect,  due  to  the  variation  of  the  current  profile  across  the  slab,  which 
leads  to  quenching  of  the  instability  if  the  slab  width  is  reduced  to  only 
a  few  ion  gyro  radii. 

In  this  paper  we  illustrate  the  importance  of  these  non-local 
treatments  of  the  CDIC1  for  the  space  plasmas.  We  show  that  the  non-local 
theory  can,  under  certain  circumstances,  lead  to  enhanced  coherence  in  the 
electrostatic  CDICI. 


II.  THEORY 

a.  Assumptions 

We  consider  a  slab  geometry  where  both  the  electrons  and  the  ions 

have  a  Gaussian  distribution  function.  The  plasma  has  a  negligible 
fiimkT 

8  (6  S  «  1),  and  the  electrons  drift  with  respect  to  the  ions  along 

the  magnetic  field  with  a  velocity  V^.  Both  species  have  a  finite  non-zero 

temperature  T  ■  m  v  2/2  k_,  a  denoting  the  species,  m  is  the  mass,  v  is 
a  a  a  B  a  a 

the  thermal  velocity  and  kg  is  the  Boltzmann's  constant.  We  neglect  the 

term  (kX  )2  «  1  where  k  is  the  wavevector  in  arbitrary  direction 
a 

while  X  is  the  Debye  length.  Since  for  most  space  plasmas  of  interest, 
a 

the  ion  and  electron  temperatures  are  of  the  same  order,  we  shall  restrict 
the  analysis  here  for  the  temperature  ratio  t  =  T^/Tg  ■  0(1).  The  ambient 
magnetic  field  is  in  the  £-2  plane  and  contains  a  shear  in  the  in¬ 
direction.  Let  us  first  consider  a  simple  model  for  a  sheared  magnetic 
field,  given  by 


B  *  Bq{z  +  (x/Lg)y},  (1) 

with  x/Lg  «  1  and  Lg  »  (30/3x)  1  where  6  »  tan  (By/Bz).  Thus  Lg  is  the 
scale  length  characterizing  the  magnetic  shear.  In  the  absence  of  shear 
the  field  configuration  is  B_  *  BQ2. 


The 


In  a  shear  free  field  the  plasma  is  unstable  to  the  CDICI  when 

h)  -  k  V.  <  0  for  id  and  k  consistent  with  the  dispersion  relation, 
la  t 

waves  at  maximum  growth  (for  Tg  *  T^)  are  characterized  by  cir  ~ 

8^  ■  eBQ/m^c,  the  ion  Larmour  frequency;  y  «  wr;  and  kp^  «  0(1)  where 


i/2 

Pi  -  (2Vi/ffli)  /8^,  is  the  ion  Larmour  radius.  Electron  Landau  and  ion- 
cyclotron  damping  reduce  the  growth  rate  depending  on  the  magnitude 
of  k  (see  Ganguli  and  Bakshi  (1982)). 


The  above  description  of  the  plasma  gets  significantly  modified  by  the 

introduction  of  a  shear  in  the  magnetic  field  as  given  in  expression  (1). 

The  magnetic  field  rotates  in  the  £-£  plane  (see  Figure  1)  as  a  function  of 

x.  If  at  x  ■  0  we  have  k  »  0  then  at  x  *  x,we  see  that  k  *  0  since  at  x, 

H  1  I!  * 

the  magnetic  field  has  rotated  by  an  angle  61“  x1/Lg.  Thus  the 
dispersive  properties  of  the  plasma  are  also  a  function  of  x.  It  is  this 
effect  which  changes  the  boundary  condition  from  plane  waves  to  outgoing 
energy  flux  condition  and  which  the  local  theory  fails  to  account  for, 
thereby  giving  incomplete  results. 

We  introduce  shear  (i)  locally,  by  replacing  k  by  k°  +  sk  x,  where 

m  ®  ®  y 

8  •  1/LS  and  (ii)  globally  by  replacing  il^  by  "  —  ".  Since  the  magnitude 
of  the  magnetic  shear  of  interest  to  us  is  quite  small,  we  neglect  the 
orbital  effects  of  shear  [Bakshi,  Bellew,  Ganguli  and  Satyanarayana  (1977), 
Bellew  (1978),  Ganguli  (1980)  and  Linsker  (1981)]  arising  mainly  out  of  the 
shear  kinematic  drift  of  the  particles  in  the  sheared  magnetic  field, 
b.  Dispersion  Relation 


The  general  dispersion  relation  for  CDICI  in  the  absence  of  shear 
is  given  by 


k2  +  £ 
n,  a 


T  (b  )  ur-k  V  uj-nn  -k  V, 

n  a  r>  .  c  1  da>  „  r  a  II  da^ 

“12 —  i1  +  qrrr  )  z  (— iirp —  JJ 


0, 


IT 


II 


(2) 


where  b  •  k2p2/2,  r  (b)  *  I  (b)e  ^  and  I_(b)  are  modified  Bessel's 
q  1  ot  n  n 

function. 

As  previously  described  we  study  the  effects  of  the  magnetic  shear  by 
replacing  k(  by  kyU,  where  u  ■  sx(  =  k|](x)/ky)  and  ikx  by  —  in  equation 
(2)  (we  have  assumed  k°  *  0).  Transforming  to  the  ion  frame  (i.e.,  setting 
Vdi  ■  0),  retaining  the  n  ■  0  term  for  the  electrons  and  n  *  0,  +  1,  +  2 


3 


terns  for  the  ions  (since  we  wish  to  study  the  first  harmonic  only)  and 
expanding  rn(b-(p|/2)82/3x2)  to  0(a2/3x2)  (here  b  «  k2p2/2),  we  obtain  a 
second  order  differential  equation  for  the  electrostatic  potential  $(x). 


{-A(x)  -|j7  +  Q:(x)}*(x)  «  0. 


(3) 


The  equation  can  be  reduced  to  (see  Ganguli  and  Bakshi  (1982)), 
{ ( p^s) 2  +  Q(u)}*(u)  ■  0, 


(4) 


where. 


Qi(u)  . 

Q(u)  - - : -  , 

A(u) 


Q2(u)  -  1  +  t  +  CirQZ(ci)  +  51r1{Z((l-p);i)  +  Z  ((l+p)^)} 


+  qr2{z((i-2p)Ci)  +  z((i+2p)Ci)} 

k  nV  j  k  V 

-Hr»’ 

A(u)  -  (ii)[Z(Ci)r;  +  {Z((i-p)ct))  +  Z((i4p)ct))rl 


+  {z((i-2p)Ci)  +  Z((l+2p)Ci)}r2]f 


and 


We  have  assumed  k2  «  k2  in  neglecting  higher  order  derivatives  in 

x  y 

the  expansion  of  rn» 

As  shown  in  Ganguli  and  Bakshi  (1982)  we  expand  the  "potential",  Q(u) 
of  the  differential  equation  (4)  around  £  *  (u-uo)  “  0  t0  °(£2)*  Here  uQ 
is  the  angle  of  propagation  for  the  maximum  growth  rate  in  the  absence  of 
shear.  We  thus  obtain, 

Q  " 

{(pjS)2  -|^2-  +  a  +  — n2}4>(n)  ■  0,  (5) 

where, 

"  ■ « * 


4 


Ml 


Q  -  Q(u  ),  Q'  -  —  Q(u)  ,  Q  "  -  -5— r  Q(u) 
xo  x  o'*  xo  du  x  '  u*u  *  xo  du2  x  'u“u 

o  o 


The  dispersion  relation  under  the  outgoing  energy  flux  boundary 
condition  similar  to  Pearlstein  and  Berk  (1969),  is 

Q(w,u  )  -  (2  £+1) (p. s)(-Q  ”/2) 1/2  +  Q'  2/2Q  (6) 

O  10.00 

Note  Q(u),u  )  *  0  yields  the  local  dispersion  relation  maximized  over  k  , 
while  the  right  hand  side  of  the  equation  (6)  arises  from  of  the  non-local 
treatment.  The  first  term  in  the  right  hand  side  is  proportional  to  the 
magnetic  shear  and  vanishes  in  the  zero  shear  limit.  The  second  term,  is 
not  explicitly  shear  dependent.  This  term  can  contribute  significantly 
even  for  infinitesimal  shear  thereby  making  the  non-local  dispersion 
relation  solutions  much  different  from  those  of  the  local  dispersion 
relation. 

c.  Finite  Slab  Effects. 

We  have  given  elsewhere  (Bakshi,  Ganguli  and  Palmadesso,  1983)  the 
details  of  the  modifications  brought  about  by  the  finite  slab  size  Lc.  A 
finite  width  current  is  introduced  by  taking  the  electron  distribution 
function  to  be  a  shifted  Maxwellian  with  a  drift  velocity  parallel  to  the 
magnetic  field  given  by 


V*)  -  V°  8(*g/Lc), 


8(0  -  e'5  ,  *g  -  x  *-Z  . 


n 


If  we  ignore  the  very  small  effects  of  order  p2/L2,  where  p  is  the 

e  c  e 

electron  gyro-radius,  the  current  profile  is  found  to  be 


5 


Such  a  current  profile  generates  a  self-consistent  shear  in  the  magnetic 
field ,  given  by 


B  (x) 

y 

B 


Lc  « 

I rl 

s  o 


where  the  shear  length  is  defined  by 

cB 


±«L 


4un  eV° 
o  d 


The  corresponding  variation  in  the  parallel  wavenumber, 


(9) 


(10) 


kjU)  -  k°  +  kyBy(x)/Bz,  (11) 

leads  to  the  variation  in  the  angle  of  propagation 

k  (x)  L  £ 

u(x)  - - - -u  +?T  /  8(5)dC*  (12) 

y  so 

Thus  the  relation  between  the  angle  u  and  the  physical  distance  x,  or  the 
scaled  physical  distance  £  is  no  longer  linear.  The  non-local  effects  are 
again  described  by  the  differential  equation  obtained  by  letting 
k^  (l/i)d/dx,  and  k(  ♦  kt(x)  as  given  by  eq.  (11)  and  Vd  -►  Vd(x)  as 
implied  by  eq.  (8), 

[pf  |^7  +  Q(u(x),  vd(x))]$  -  0  (13a) 


or 

KirO2  +  ■  °»  (13b> 

c  ’ 

with  Q  defined  as  in  (4)  apart  from  the  x-dependence  as  stated  above. 
Expanding  Q  around  where 


provides  more  accurate  results  as  compared  to  expansion  around  or  uQ  as 
in  eq.  (5),  and  the  dispersion  relation  which  generalizes  eq.  (6)  is  given 
by 

QU,^)  -  (2 1  +  l)(Pi/Lc)[-  Y  (TU,^)]172.  (15) 

The  eigenf requency  obtained  from  (15)  depends  on  k°,  which  is  as  yet 
undetermined.  The  dominant  mode,  for  a  given  ky,  is  obtained  by 

maximizing  y  *  Ima)(k .)  over  k°. 

II  II 


III.  RESULTS 

We  now  discuss  the  implications  of  the  non-local  treatments  in  the 
preceding  section.  The  results  of  the  uniform  shear  model  are  contained  in 
eq.  (6),  while  the  finite  slab  effects  are  given  by  eq.  (15).  Both  these 
are  generalizations  of  the  local  dispersion  relation  Q(u,uq)  *  0. 

The  non-local  dispersion  relation  (6)  is  solved  numerically  by  an 
iterative  technique.  We  simplify  the  "potential"  by  expanding  all  the  ion 
Z-f unctions  except  n«l  for  large  arguments,  retain  terms  to  0(l/^2)  and 
treat  A(u)  as  a  constant  evaluated  at  uQ.  [For  details  see  Ganguli  and 
Bakshi  (1982)]. 

Figure  2  shows  a  plot  of  y/Q  ^  against  b  for  various  values  of  the 
shear  and  r  *  0.5,  y  ■  1837,  V^/ve  -  0.25.  Note  the  difference  between  the 
local  solution  and  the  zero  shear  limit  of  the  non-local  solution.  This 
large  reduction  in  the  growth  rate  is  due  to  the  second  term  in  equation 
(6)  arising  from  the  non-local  treatment.  Figure  2  also  shows  that  the 
value  of  b  for  which  maximum  growth  occurs  remains  essentially  unaltered 
from  the  local  theory,  even  when  shear  is  increased.  In  addition  to  the 
expected  reduction  of  the  maximum  growth  rate,  the  primary  effect  of 
Increased  shear  is  to  narrow  the  band  of  unstable  perpendicular 
wavelengths.  The  uneven  spacing  of  the  p^s  *  0  ,  0.005  and  0.01  curves 
shows  that  non-linearity  in  the  shear  parameter  (through  the  implicit 
dependence  of  u  on  (p1s))already  sets  in  at  o^s  -  0.01. 


7 


In  Figures  3  and  4  we  provide  additional  y/n  ^  versus  b  plots  for 

typical  space  plasma  parameters.  In  Figure  3  we  consider  a  hydrogen 

_6 

plasma:  p  *  1837,  t  ■  0.75,  V^/Vg  "  0*3  and  p^s  *  10  .  Note  once  again 

the  big  drop  in  the  local  growth  rate  as  well  as  the  significant  reduction 

of  the  unstable  perpendicular  wavevector  band.  Figure  4  illustrates  the 
same  features  for  t  **  1,  V^/Vg  =  0.35  and  the  other  parameters  unchanged. 
The  consequence  of  the  non-local  treatment  is  even  more  drastic  in  this 
case.  The  non-local  treatment  leads  to  a  completely  damped  instability. 

The  drastic  changes  due  to  non-local  effects  stem  from  the  implicit 

assumption  that  the  driving  current  prevails  over  a  sufficiently  wide 

region  in  the  x-direction  so  that  k  /k  can  become  arbitrarily  large.  This 

n  -L 

condition  will  not  be  met  in  many  space  plasma  situations,  and  we  need  the 
results  of  the  finite  slab  theory  (Bakshi,  Ganguli,  and  Palmadesso,  1983), 
which  led  to  eq.  (15). 

Consideration  of  a  finite  width  current  sheet  leads  to  two  scale 

lengths,  Lg  and  Lc,  which  govern  the  growth  rate.  Figure  5  shows  the 

variation  in  (y/j^)  vs.  (l^/p^,)  for  t  *  0.5,  V°/v  *  0.28,  p  *  1837  and  b  * 

0.6,  for  three  different  shear  parameters  (p. /L  )  *  10“*,  10“5  and  10“6. 

1  s 

In  each  case  the  growth  rate  has  the  significantly  reduced  value, 
characteristic  of  the  uniform  shear  model,  if  Lc  is  taken  to  be 
sufficiently  large.  As  the  slab  width  is  decreased,  the  growth  rate 
increases  and  climbs  up  to  the  local  value  when  Lc  is  sufficiently  small. 
Thus,  there  is  a  smooth  transition  from  the  pure  shear  reduced  y  to  the 
usual  local  y. 

The  three  curves  in  Figure  5  appear  very  similar,  except  for  a 
relative  displacement,  and  the  essential  universality  of  the  transition 
curve  (for  small  shears)  can  be  brought  out  by  plotting  y  vs.  (L_/L_)  as  in 
Figure  6.  The  shear  reduced  y  prevails  for  L„  >  L_,  and  the  three  curves 

w  v  9 

continue  to  overlap  for  Lc  >  (Lg/20).  The  curve  describing  the  weakest 

shear  p, /L  *  10  ^  almost  attains  the  local  value  when  L„  »  10“^  L_.  The 
is  c  s 

physical  reason  for  this  transition  to  the  local  value  is  easily  seen:  the 
detuning  or  misalignment  of  the  magnetic  field  lines  across  the  length  of 
the  slab  of  order  Lc  here  happens  to  be  too  small  to  have  any  impact  on  the 
instability.  On  the  other  hand,  when  Lc  >  Lg,  enough  misalignment  occurs 
to  produce  local  regions  where  the  angle  between  the  fixed  wave  direction 


8 


it  and  the  varying  S(x)  is  sufficiently  different  fron  the  optimal  value 

(typically  k°/k  -  1/10)  to  produce  a  much  lower  average  growth  rate. 

11^ 

An  additional  non-local  effect,  directly  due  to  the  variable  current 
profile  (and  not  the  shear)  is  also  evident  in  Figures  5  and  6.  After 
attaining  the  local  growth  rate,  if  the  slab  size  is  further  reduced,  the 
growth  rate  again  becomes  smaller  than  the  local  value  for  <  102. 

It  is  clear  from  Figure  5  that  this  effect  is  not  sensitive  to  the  value  of 
the  shear  (at  least  in  the  range  considered).  The  physical  explanation  for 
this  reduction  is  also  easily  seen.  The  spatial  variation  of  the  current 
has  a  direct  effect  on  the  wave  packet  The  reduction  of  the  current  away 
from  the  center  of  the  slab  results  in  a  reduced  drive  and  a  corresponding 
reduced  y.If  Lc/ pi  is  further  reduced  to  10  or  lower,  a  significant 
reduction  y  occurs,  with  eventual  quenching. 


A  typical  potential  function  4)  is  displayed  in  Figure  7,  with  V°/vg  * 
0.28,  t  -  0.5,  v  -  1837,  b  -  0.6,  L  /Pi  -  102  and  L  /Pi  «  ».  The 

independent  variable  u  is  the  physical  angle  u  »  x/Lg,  the  origin  of  x  is 
at  the  position  where  the  field  line  is  perpendicular  to  the  given  £.  As 
shown,  |<|>|  is  significant  between  0.06  and  0.12  with  a  peak  around  umo  * 
0.09.  This  indicates  that  the  wave  propagates  in  the  direction  given 

by  k°  -  (0.09)k  . 

8  X 

Besides  the  changes  in  growth  rate,  the  non-local  theory  also  predicts 
a  transition  in  the  angular  position  of  the  peak  of  |  <p  | ,  or  in  the 
effective  angle  umQ  between  the  direction  normal  to  the  local  field  line 
and  the  direction  of  the  propagating  wave.  This  is  displayed  in  Figure  8, 
for  the  same  parameters  as  in  Figure  6,  over  a  wide  range  of  (Lc/Lg).  Note 
the  transition  from  u^  ■  0.13  characteristic  of  local  theory  to  the 
reduced  *  0.09  of  the  pure  shear  dominated  limit.  The  non-local 

effects  thus  reduce  the  effective  angle  k°/k,. 

n  1 

Figure  9  shows  the  variation  of  y  as  a  function  of  wavenumber.  The 
local  theory  result  prevails  for  weak  shear  and  L„/L.  <  10  ,  the  pure 

shear  dominated  reduction  occurs  for  L„  >  L_  with  a  smooth  transition  in 
between.  The  Lc  ■  Lg  curves,  except  for  a  slight  change  of  parameters, 
correspond  to  the  p^s  *  0  and  p^s  ■  0.01  curves  of  Figure  2. 

We  have  also  examined  the  CDICI  in  an  oxygen  plasma  and  note  that  the 
current  thresholds  are  generally  much  lower  compared  with  those  of  a 


V 
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hydrogen  plasma.  The  effect  of  the  non-local  treatment  for  a  multispecies 
plasma  is  an  interesting  question  which  is  now  being  investigated  and  will 
be  reported  elsewhere. 

IV.  DISCUSSION 

There  are  several  implications  of  the  non-local  effects  discussed 
above.  If  the  shear  is  small,  we  have  shown  that  there  are  three  different 
regimes  depending  on  the  scale  size  of  the  current  channel,  Lc: 

(1)  For  Lc  >  (Ls/10)  there  occurs,  as  compared  to  the  results  of  the 
local  theory,  (a)  a  reduction  in  the  growth  rate,  (b)  a  reduction  in 
the  bandwidth  of  over  which  the  instability  prevails  and  (c)  a 
reduction  in  the  angle  of  propagation  umQ  »  k°/k^.  This  is  the  shear 
dominant  regime. 

(2)  For  p^  «  «  Lg,  one  recovers  the  results  of  the  local  theory 

(except  for  one  important  aspect,  which  we  discuss  later  in  this 
section). 

(3)  For  Lc  >  pi,  again  (a)  there  is  a  reduction  in  the  growth  rate, 
and  (b)  a  reduction  in  the  bandwidth  of  k^  over  which  the  instability 
prevails,  but  (c)  there  is  no  reduction  in  u^,  as  compared  to  the 
results  of  the  local  theory.  The  instability  is  quenched  for 

Lc  <  2er 

Thus  we  have  shown  that  local  theory  cannot  be  applied  without 
reservations;  one  must  first  scrutinize  the  scale  lengths  for  the  given 
laboratory  experiment,  or  the  space  plasma  of  interest,  to  ascertain  which 
regime  is  relevant.  In  some  situations,  (e.g.,  the  case  Ls  -  10  2p,  of 
Figure  6),  the  conditions  for  regime  (2)  cannot  be  strictly  satisfied  for 
any  value  of  Lc,  and  the  results  of  the  local  theory  are  not  applicable. 
This  aspect  is  accentuated  as  the  shear  Increases. 

In  space  plasmas,  typical  shear  lengths  are  Lg  ~  few  hundred  km,  and 
p^/Lg  »  10  Thus  the  smaller  sized  (-1  km)  auroral  arcs  would  not 
qualify  for  regime  (1).  The  quiet  time,  broader  Birkeland  current 
Region  1/Region  2  systems  may  be  more  appropriate  for  these  effects.  Also, 
since  the  current  sheet  width  increses  with  altitude,  the  effects  of  regime 
(1)  have  a  better  chance  of  being  significant  at  higher  altitudes. 


The  transition  from  the  shear-affected  regime  (1)  to  the  local  theory 

results  of  regime  (2)  was  shown  (Figure  6)  to  depend  on  L  /L_,  essentially 

c  s 

described  by  a  universal  curve  independent  of  shear.  The  net  change  in  he 
growth  rate  and  the  range  of  (Lc/Lg)  over  which  the  transition  takes  place 
will  be  governed  by  the  physical  parameters  p,  t  and  V^/ve,  and  the 
perpendicular  wavenumber  k^.  Under  what  conditions  can  we  expect  the 
effects  of  shear  to  be  significant? 

One  can  assess  these  trends  by  plotting  the  local  growth  rate  as  a 
function  of  the  angle  k^/k^.  As  an  illustration,  we  note  that  Figure  1  of 
Ganguli  and  Bakshi  (1982)  displays  the  effect  of  varying  t.  It  shows  that 
the  angular  band  over  which  the  instability  can  occur  becomes  narrower  as 
t  increases.  The  transition  from  regime  (1)  to  regime  (2)  can  be  expected 
to  occur  when  Lc/Lg,  which  represents  the  tilt  produced  by  shear,  becomes 
comparable  to  this  angular  bandwidth.  Thus  one  can  infer  that  for  a  given 
Lc,  increasing  x  will  favor  regime  (1),  making  the  shear  more  effective. 

We  can  test  this  quantitatively  by  evaluating  the  growth  rate  as  a 
function  of  L  /L_  for  two  different  values  of  x.  Figure  1  of  Ganguli  and 
Bakshi  (1982)  was  for  a  fixed  b  *  0.5.  Since  the  perpendicular  wavenumber 
(or  equivalently  b)  is  a  free  parameter,  one  must  maximize  the  growth  rate 
over  b  as  well.  We  have  carried  out  these  calculations  for  V^/Vg  •  0.4, 
p  «  1837,  t  ■  0.5  and  1,  and  the  results  are  displayed  in  Figures  10  and 
11.  The  imaginary  part  of  Qj(=  Q^),  e9*  (4),  is  plotted  against  u  * 
k|/kj.  *’*8ure  with  b  -  0.6,  coo  *  1.3037  for  t  ■  0.5  and  b  *  0.9, 

u>0  ■  1.2027  for  t  ■  1.  The  angular  width  of  the  domain  where  Qj^  is 
negative,  indicative  of  instability,  is  clearly  much  smaller  for  x  *  1. 

The  corresponding  results  for  the  growth  rate  as  a  function  of  (Lc/Lg) 
are  displayed  in  Figure  11.  The  growth  rate  for  each  x  is  normalized  to 
the  corresponding  local  growth  rate  yT  •  Noticeable  departure  from  the 
local  value  begins  at  (L_/L_)  *  0.01  for  x  *  1  and  at  0.02  for  x  *  0.5,  and 
proceeds  much  faster  for  the  larger  x  as  Lc/Lg  is  increased.  Thus  the 
effects  of  shear  are  more  significant  for  the  larger  x,  given  a  fixed 
current  channel  width. 


Figure  1  of  Ganguli  and  Bakshi  (1982)  also  shows  that  the  angular  band 
width  of  the  domain  of  instability  increases  with  increasing  V^/Vg.  But  an 
increase  in  also  strengthens  the  shear,  thereby  increasing  Lc/Ls.  Near 
the  threshold  drift  velocity,  (where  the  system  is  marginally  stable),  the 
increase  in  angular  band  width  proceeds  faster  than  the  increase  in  Lc/Lg; 
thus  an  increase  in  favors  a  move  towards  the  local  regime  (2).  On  the 
other  hand,  if  is  significantly  larger  than  the  threshold  for  marginal 
stability,  the  situation  is  reversed  and  an  increase  in  favors  a  move 
towards  the  shear  dominated  regime  (1). 

So  far,  we  have  only  discussed  the  mode  near  the  cyclotron  frequency, 
u)  >  u>ci.  What  can  we  expect  for  the  higher  harmonic  modes?  The  definition 
of  Qj  in  eq.  (4)  must  first  be  modified  to  include  higher  n  terms  from  the 
full  dispersion  relation,  eq.  (2).  We  have  evaluated  Qj  including  terms  up 
to  n  *  +  4  and  maximized  over  b  to  obtain  optimal  local  growth  rates.  The 
results  for  Qjj  vs.  u  are  displayed  in  Figure  12  for  the  first  three 
cyclotron  harmonics,  with  real  frequencies  u>  ■  1.1897,  b  ■  0.9;  u)q  * 
2.2237,  b  *  2.4;  u  -  3.2499,  b  «  4.8;  and  physical  parameters 
t  *  1,  y  *  1837,  V^/Vg  -  0.65.  The  angular  band  width  of  the  unstable  zone 
(Qll  <  0)  becomes  significantly  smaller  for  higher  harmonics.  Thus  the 
transition  to  the  shear  dominated  regime  (1)  can  be  expected  to  occur  for 
the  higher  harmonics  for  progressively  smaller  values  of  Lc/Lg. 

These  considerations  show  that  the  results  of  the  local  theory,  regime 
(2),  cannot  be  used  without  reservations,  and  a  transition  towards  regime 
(1)  is  helped  by  increasing  t,  considering  higher  harmonics,  and  increasing 
the  drift  velocities.  It  is  also  accentuated  by  increasing  the  current 
channel  width,  directly  by  considering  wider  channels  or  for  the  same 
current  sheet,  by  scaling  with  altitude.  A  detailed  numerical  study  of  all 
of  these  effects  over  a  wide  range  of  parameters  is  beyond  the  scope  of 
this  work  and  will  be  given  elsewhere.  We  anticipate  that  some  of  the 
situations  of  practical  interest  in  space  plasmas  will  fall  in  regime  (2) 
while  others  will  favor  a  move  towards  regime  (1). 

There  is  a  subtler  effect,  implicit  in  the  non-local  treatment,  which 
may  be  more  widely  relevant.  Once  any  inhomogeneity  is  introduced  in  the 
x-direction,  the  symmetry  in  the  x-y  plane  is  broken,  and  one  has  to 

A 

distinguish  between  the  propagation  along  y  and  the  propagation  (if  any) 


along  x.  The  same  inhomogeneity  precludes  the^se  of  Fourier  transforms  in 

x 

the  x-direction  and  the  plane  wave  solutions  e  are  necessarily  replaced 
by  a  wave  packet,  $(x),  which  is  localized  (e.g.,  see  Fig.  7)  in  the  x- 
directlon.  Since  the  current  systems  have  the  structure  of  a  sheet,  the 


short  direction  becomes  the  slab  direction  and  the  waves  propagate  along  y 

A 

and  z  with  wave  vector  components  (ky,  V-  ky  »  V 


Now  let  us  consider  the  implications  of  this  phenomenon  concerning 
satellite  measurements  of  ion  cyclotron  waves .  A  moving  observer  will  see 
a  Doppler  shifted  frequency 


(o'  ■  a)  +  k.v. 


where  v^  is  the  velocity  vector  of  the  observer  (satellite)  and  _k  is  the 
wave  vector.  Since  k_  has  no  component  along  x  (the  wave  packet  is 
localized),  it  follows  that  a  satellite  moving  perpendicular  to  the  current 
sheet  (i.e.,  along  our  x  direction)  will  not  see  any  Doppler-shif t, 
irrespective  of  the  size  of  ky.  Such  a  statement  could  not  be  made  in  the 
purely  local  theory,  since  there  would  be  no  essential  difference  between 

A  * 

the  two  transverse  directions  x  and  y,  and  waves  are  as  likely  to  be  moving 
in  any  direction  in  the  x-y  plane.  In  a  local  theory,  a  wide  range  of 
Doppler  shifts  would  be  produced,  from  -k^vg  to  +k^vg,  since  the  maximum 
value  of  is  k^.  Some  of  the  observations  on  S3-3  indicate  no  or  little 
Doppler  shift*  While  other  explanations  may  be  available,  this  inference 
from  our  non-local  treatment  must  also  be  considered,  since  an 
inhomogeneous  plasma  mandates  a  non-local  treatment.  Thus,  the  extremely 
high  degree  of  coherence  of  many  of  the  observations  can  be  understood  in  a 
very  simple  way.  This  idea  can  be  further  tested  by  analyzing  satellite 
data  from  an  orbit  which  crosses  the  current  sheet  at  a  significantly 
different  angle. 

Various  phenomena  described  above  would  lead  to  an  enhanced  coherence 
for  the  observed  ion  cyclotron  waves  as  seen  by  a  moving  observer.  When 
conditions  for  regime  (1)  prevail,  there  is  a  significant  reduction  in  the 
band  of  perpendicular  wave  numbers  for  which  the  instability  can  occur. 
This  reduces  the  magnitude  of  the  Doppler  shift  k.v^  as  compared  to  the 
local  theory  for  any  direction  of  the  satellite  motion.  On  the  other  hand, 
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the  directional  preference  for  It  ,  required  by  the  inhomogeneity  for  all 
regimes,  shows  that  satellites  moving  perpendicular  to  the  current  sheet 
would  see  very  coherent  ion  cyclotron  waves,  with  no  Doppler  shift, 
irrespective  of  the  size  of  A  study  of  all  available  observational 
information  for  various  transit  angles  of  the  satellite  motion  could  be 
very  Instructive. 

V.  OUTLOOK 

In  this  paper,  we  have  analyzed  the  implications  of  a  non-local 
treatment  of  the  CDIC1,  which  is  necessitated  by  the  introduction  of  an 
inhomogeneous  current  profile.  One  cannot  yet  provide  a  detailed 
comparison  between  theory  and  observations;  that  would  require  additional 
work  on  the  theoretical  side  and  more  detailed  data  from  the  observational 
side.  New  data  from  the  Dynamic  Explorer  would  be  very  helpful  in  this 
context . 

But  even  without  carrying  out  the  additional  theoretical  studies,  or 
awaiting  the  new  observational  data,  a  number  of  interesting  questions  have 
been  raised  in  this  paper,  and  it  would  be  very  helpful  to  look  at  all 
available  data  from  a  different  angle,  in  terms  of  the  framework  of  this 
paper. 

The  model  of  inhomogeneity  of  the  current  employed  here  was  one¬ 
dimensional.  Further  generalizations  of  theory  are  needed  if  the  current 
sheets  are  of  limited  extent  in  both  transverse  directions.  Even  within  a 
one-dimensional  model,  there  could  be  a  sub-structure  on  a  shorter  distance 
scale,  or  even  a  random  or  stochastic  current  structure  superposed  on  the 
underlying  smooth  current  profile.  These  aspects  also  would  require 
further  generalizations.  Thus  it  would  be  very  valuable  to  learn  from  the 
present  and  future  observations  the  detailed  nature  of  the  current 
distributions. 

It  is  equally  important  to  classify  all  the  available  observations  on 
the  ion  cyclotron  waves  according  to  the  satellite  transit  direction  and 
study  its  correlations,  if  any,  to  the  peak  frequency  and  the  width  of  the 
observed  power  spectrum. 
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Further  improvements  in  understanding  the  CDICI  can  come  in  a  number 
of  other  ways  as  well.  Additional  generalizations  of  theory  should  include 
variation  (i.e.,  inhomogeneities)  along  the  field  lines.  In  terms  of 
applications  of  the  theory  of  the  present  level  one  should  certainly 
consider  multi-species  plasmas,  wider  ranges  of  the  ion  to  electron 
temperature  ratio  t,  wider  ranges  of  the  driving  current  V°/vg,  and  effects 
on  higher  harmonics.  In  terms  of  interpretation  the  results  can  be  recast 
in  terms  of  the  convective  growth  rate  approach  (Ashour-Abdalla  and  Thorne, 
1978),  which  takes  into  account  the  small  group  velocity  of  the  ion 
cyclotron  waves.  These  and  other  related  topics  are  under  further 
investigation. 

We  also  note  that  Hwang,  Fontheim  and  Ong  (1983)  have  recently 
considered  the  effects  on  electrostatic  waves  of  a  cold  electron  current 
sheet  of  finite  thickness,  and  concluded  that  the  current  sheet  partially 
stabilizes  the  system  and  contributes  to  the  coherence  of  the  excited 
waves.  Since  they  did  not  consider  the  effects  of  magnetic  shear,  the 
correspondence  is  with  our  regime  (3),  where  the  non-local  effects  stem 
directly  from  the  finite  width  of  the  current  channel.  Even  in  this  regime 
their  approach  is  restricted  to  «  L£.  Also  their  assumption  of  a  sharp 
change  in  the  current  at  the  boundaries  is  not  pertinent  to  the  space 
plasmas.  We  also  note  that  they  have  concentrated  more  on  the  ion-acoustic 
wave,  and  have  no  thermal  spreads,  but  some  of  their  conclusions  are 
similar  to  ours. 
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Figure  1.  A  sketch  of  a  sheared  magnetic  field.  At  x  -  0  the  magnetic 
field  BQ  is  in  the  i-direction.  At  x  ■  Xj  the  magnetic  field 
rotates  by  an  angle  0^  ■  sx^,  where  s  *  1/LS.  A  wave 

perpendicular  to  the  magnetic  field  B0  at  x  ■  0  is  no  longer  so 


at  x  «  Xi 


Figure  2.  A  plot  of  the  growth  rate  (y/f^)  against  b  for  various  shear 
values.  Here  t  ■  0.5,  y  ■  1837,  and  V^/ve  *  0.25.  The  big 
difference  between  the  local  solution  and  the  (p^s)  ■  0  limit 
of  the  non-local  solution  is  due  to  a  shear  independent  term 
arising  from  the  non-local  treatment.  The  curves  also  indicate 
that  the  effect  of  the  shear  is  more  pronounced  for  larger  b. 


NONLOCAL 


A  plot  similiar  to  Figure  2.  Here  t  *  0.75,  y  -  1837  and 
V^/vg  ■  0.3.  We  provide  only  the  local  and  the  non-local 
spectrum.  The  non-local  curve  has  a  very  small  shear, 

(p^s)  ■  10  Note  the  significant  reduction  in  the  unstable 
perpendicular  wavevector  band. 


Figure  4.  Growth  rate  versus  b  plot  for  y  *  1837,  and  V^/Vg  ■  0.35.  The 
effect  of  the  non-local  treatment  is  very  drastic  in  this 
case.  The  instability  is  totally  suppressed. 
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Figure  5.  A  plot  of  the  growth  rate  (•^-)  against  the  current  channel 
width  (kg/PjJ  for  three  different  magnetic  shear  lengths.  Here 
V°  -  0.28  v  .  u  ■  1837,  b  «*  0.6,  and  t  *  0.5. 


Figure  6.  A  plot  of  against  (Lc/Lg)  for  three  different  Lg 

values.  Here  b  ■  0.6,  y  ■  1837,  t  *  0.5,  and  V°  ■  0.28  ve> 

The  transition  from  local  to  the  nonlocal  results  are 
universal,  almost  independent  of  the  shear  value. 
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Figure  8.  A  plot  for  the  real  angle  Umo  where  |$(u)|  attains  its  maximum, 
against  Lc/Lg.  The  physical  parameters  are  same  as  in 
Figure  6. 


Figure  9.  A  plot  of  the  growth  rates  against  b.  Here  t  ■  0.5,  p  *  1837, 

V°  «  0.28  v  .  The  dotted  lines  show  the  transition  from  local 
a  e 

to  nonlocal  values  as  L„  approaches  L_.  The  lowest  solid  line 

2 

is  for  L.  *  L_  ■  10  o . . 


Figure  10.  A  plot  of  vs.  u  for  t  *  0.5  and  1,  maximized  over  b  to 

attain  maximal  growth  rates. 
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Figure  11.  A  plot  of  y/yl  vs.  Lc/Lg  for  t 
of  Figure  10. 


0.5  and  1,  for  the  parameters 


Figure  12.  A  plot  of  Qjj  vs.  u  for  the  firBt  three  cyclotron  harmonics  for 
r  •  1,  Vjj/Vg  ■  0.65,  maximized  over  b  to  attain  maximal  growth 
rates.  The  solid  line  is  the  third  harmonic,  the  dashed  line 
with  small  dashes  is  the  second  harmonic  and  the  dashed  line 
with  the  large  dashes  is  the  first  harmonic. 
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